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Abstract 

We prove that a gauge approach based on a composite struc¬ 
ture of the hole in hole-doped cuprates is able to capture analyti¬ 
cally many features of the experimental data on superfluid density in 
the moderate-underdoping to nearly-optimal doping region, including 
critical exponent, the Uemura relation and near universality of the 
normalized superfluid density. 


1 Introduction 

The superfluid density ps of hole-doped cuprates in the region of moderate 
underdoping to nearly optimal doping as a function of temperature exhibits 
features whose combination does not allow for a simple explanation. The 
well-defined gapped Fermi arcs seen in ARPES experiments [1] and the linear 
T-dependence near T = 0 of hints at a BGS-like hole dynamics. Attribut¬ 
ing the low temperature T-linear dependence to the quasi-particle excitations 
near the nodes of the d-wave BGS order parameter [2], however, it is diffi¬ 
cult to reconcile with the non-mean-field critical exponent 2/3, typical of 
a 3DXY model, at the superconducting transition [3j(in a BGS theory one 
would expect the mean-field value 1). Also it appears at odds with the in¬ 
triguing Uemura linear relation |1] between p^iT = 0) and approximately 
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verified for moderate underdoping (in a BCS theory the superfluid density at 
T = 0 does not depend on the BCS order parameter controlling T^). An alter¬ 
native explanation of the T-linear behaviour is based on phase [5] or pairing 
fluctuations |6] of the order parameter in a BCS-BEC crossover setting. In 
both cases one can obtain a low-energy effective theory in terms of the phase 
field described by an XY model, whose superfluid density exhibits the desired 
linear slope. Typically the advocated XY model is two-dimensional, but then 
the critical behaviour comes out experimentally incorrect. A 3-dimensional 
nature of the XY model is sometimes claimed to emerge in a narrow range 
of temperatures close to Tc due to the presence of a stack of Cu-0 layers 
in the high-T^ materials. However the range of a 3-dimensional behaviour 
found experimentally is larger than expected and when one crosses to the 
region where the 2-dimensional behaviour is anticipated one does not find 
the jump at the Kosterlitz-Nelson temperature, characteristic of superfluid 
density in the BKT formalism[7]. Moreover, it is possible to interpret the 
large pseudogap experimentally observed above in terms of the BCS-BEC 
crossover, however it would be at odds with the presence of a well defined 
gapped Fermi surface. 

In this paper we show that a solution to the above puzzles is provided by 
a gauge approach developed in [8] [9] [10], based on a modelling of the low- 
energy physics of Cu-0 planes by a planar t — t' — J model with square lattice 
and heavily relying on a composite nature of the hole. Such model is widely 
believed to be appropriate for the description of hole-doped cuprates due to 
the large charge transfer gap in these materials, with Zhang-Rice singlets as 
charge carriers, see e.g. m- The hole in our approach to this model is seen 
as a composite of a spin (spinon) and a charge (holon) excitation. As dis¬ 
cussed below the specific choice adopted for these excitations is different from 
the standard slave-boson approach and it is suggested by ID results. Specif¬ 
ically the charge degrees of freedom are mainly responsible for the Fermi 
surface structure, whereas in the considered doping region the superfluid 
density is dominated by the spin excitations which are also controlling the 
superconducting transition. Since is not determined by the charge degrees 
of freedom, this dichotomy allows for the existence of a well defined Fermi 
surface along with both a large pseudogap region, needed for a pairing tem¬ 
perature for the charges well above the superconducting one, and with the 
strongly BEC behaviour, needed for the XY-like superfluid density and the 
Uemura relation. Furthermore, the spin excitation dynamics is “relativistic” 
with space and time treated on the same footing and it is described by a 3D 
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XY model with temperature dependent coupling constant, thus explaining 
the critical exponent of the superfluid density. Actually it turns out that, 
consistently with the experiments in the moderate underdoping to nearly op¬ 
timal doping region, the superfluid density normalized by its value at T = 0 
as a function of T/Tc is well fitted by the 3DXY curve almost universally. 
Furthermore an approximate Uemura relation is analytically derived. 


2 The gauge approach 

We start sketching the key ideas of our gauge approach relevant to understand 
the treatment of the superfluid density. The no-double occupation of the the 
t — t' — J model is tackled with a spin-charge decomposition of the hole c 
written as a product of a spinon s, spin 1/2 boson, times a holon h, spinless 
fermion. Since h is spinless, by Pauli principle the no-double occupation 
constraint is automatically implemented. This decomposition introduces an 
unphysical degree of freedom due to a local U{1) gauge invariance, because 
one can multiply at each site the spinon and the holon by arbitrary opposite 
phase factors leaving the physical hole field unchanged. This invariance is 
made manifest with the introduction of a slave-particle gauge held A, which 
in turn produces an attraction between s and h. A long-wavelength treatment 
of the spinons leads to a CP^ spinon nonlinear a—model with A as gauge 
held, producing a “relativistic” dispersion relation for the spinons. In our 
approach we further use a feature specihc to 2 (or 1) space dimensions: we 
add a 1/2 charge-hux to h and a 1/2 spin-hux to s, still retaining the 
fermionic statistics of c [I2]. The spinless holons dressed by are expected 
to obey Haldane-Wu semion statistics [13] allowing double occupation in 
momentum-space, as can be proved for their 1-dimensional analogue, where 
the statistics is crucial to get the correct critical exponents [TT|. We assume 
this property for our dressed holons [T5| and therefore these holons have 
a Fermi surface at low T coinciding with that of spin 1/2 fermions of the 
same density. The physical hole is obtained as a holon-spinon bound state 
produced by the gauge attraction and it inherits the features of the holon 
Fermi surface, but renormalized and strongly overdamped by the interaction 
of the spinon with gauge fluctuations. [8] 

We consider a mean-held (MF) approach where we neglect the holon 
huctuation in and the spinon huctuations in The result is for the 
spin-hux ^sU) = (<^^/2) I]iarg(j - i)/i*/ii(-l)l*l , where j,i are sites of the 
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lattice and |i| = ix + iy ■ The gradient V<hs can be viewed as the potential of 
U{1) spin-vortices centered at holon positions with (—1)'*' chirality . Hence 
the empty sites of the model where holons sit are surrounded by spin-vortices, 
a quantum distortion of the antiferromagnetic (AF) spin background, with 
opposite chirality in the two Neel sublattices. These vortices are missed in 
the slave-boson treatment. The interaction of spinons and spin-vortices in 
the long wavelength limit is of the form 

J(1-25)(V$"(a;))^s*s(a;) (1) 

and it is the source of both short-range AF and charge pairing. From a 
quenched treatment of spin-vortices we derive the MF expectation value 
((V<h®(a;))^) = ml ^ 0.5(5| log(5|, which opens a mass gap for the spinon, 
consistent with AF correlation length at small 5 extracted from the neu¬ 
tron experiments [IB]. Thus, propagating in the gas of slowly moving spin- 
vortices, the AF spinons, originally gapless in the undoped Heisenberg model, 
acquire a hnite gap, leading to a short range AF order. By averaging instead 
the spinons in Eq. ([1]), we obtain an effective interaction: 

J(1 - 2i)(s-s) - mhji-h,, (2) 

id 

where A is the 2D lattice Laplacian. From Eq. ([2|) we see the interaction 
mediated by spin-vortices on holons is of 2D Coulomb type. From the known 
behaviour of planar Coulomb systems we derive that below a crossover tem¬ 
perature Tph J(1 — 25){s*s) (identihable with T* in [T7] i a hnite density 
of incoherent holon pairs appears. Therefore the origin of the charge-pairing 
is magnetic since it involves the Neel chirality structure of vortices , but it 
is not due to exchange of AF spin huctuations. In MF for the charge-hux 
we get <Fh(j) = (1/2) ~ *) yielding a vr hux per plaquette converting 

the low-energy modes of the spinless holons h into Dirac fermions with linear 
dispersion and two small Fermi surfaces (FS) centered at (±7r/2,7r/2) in the 
the magnetic Brillouin Zone (BZ) with Fermi energy ep ^ t6, characterizing 
what we name the “pseudogap phase” (PC) of the model. By increasing dop¬ 
ing or temperature one reaches a crossover line T* < Tph, identihed with the 
experimental inhection point of in-plane resistivity (and with T** in [T7]). 
then entering the “strange metal phase” (SM) in which the effect of the 
charge flux is screened by spinon phases and we recover a “large” FS for the 
holons with eiT’ ~ t(l -|- 5). In PC in the BCS approximation 2 p-wave orders 
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for the two holon Fermi surfaces recombine to give a d- wave order in the 
full BZ for the hole|18j.jl^j- In SM but below Tph the situation is slightly 
more involved, but qualitatively similar HD]. The pairing turns out to be 
not far from the BCS-BEC crossover, but still in the BCS side with kp^o 
ranging from to 1.5 to 3.5 for the doping region we considered, being the 
BCS coherence length. However, we do not have yet condensation of holon 
pairs because the fluctuations of the phase of the pairing held destroy the 
holon-pair coherence of BCS approximation and produces the phenomenol¬ 
ogy of Fermi arcs coexisting with gap in the antinodal region [10]. Charge 
pairing alone is not yet hole-pairing, since the spins are still unpaired. Only 
at a lower temperature Tpg using the holon-pairs as sources of effective spin 
attraction, the gauge attraction between spin and charge induces the forma¬ 
tion of incoherent RVB spin pairs. This attraction is energetically favorable 
since it lowers the spinon gap from to M k, — (As)^/ms , where 
is the modulus of the spinon-pair held (assumed spatially constant in MF). 
Below Tps the low-energy ehective action obtained integrating out spinons is 
a Maxwell-gauged 3D XY model, where the angle-held (j) of the XY model 
is the “square root” of the phase of the long-wave limit of the hole-pair held 

and the gauge held is A. The holon contribution is 
QED-like but subleading. Introducing integer-valued vector currents rip that 
allow (j) to be treated self-consistently as an angle function with range in the 
interval from 0 to 271 (lattice regularization is understood in the following), 
the low-energy ehective action S'efr is given [12] by: 



(3) 


where H^^^ is the vacuum polarization of the holons in the low-energy 
limit. Since in the range considered Vp/T > Vg/T ~ J/T ^ 1, with vp the 
Fermi velocity of holons and Vg the spinon velocity, we extended the integra¬ 
tion to the full Euclidean space-time retaining though the temperature 
dependence of the coefficients in the action. Neglecting the subleading con¬ 
tribution of holons the spinon pairing parameter A^ satishes a gap equation 
of the form 
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where the energy is measnred in nnits of J, E± = \JkP‘ + ml ± 2fcAs and 
r ~ |A/j|, with the holon pairing parameter, derived in the BCS approx¬ 
imation. The ganged XY model has two phases: Conlomb and Higgs. If in 
(]5j) the coefhcient of the Anderson-Higgs mass term for the gange held A, 
~ A^, is snfhciently small, the phase held (p hnctnates so strongly that no 
mass gap for A can be prodnced. This is the Coulomb phase, where a plasma 
of magnetic vortices-antivortices appears. For a sufficiently large coefficient 
of the mass the gauged XY is in the broken-symmetry phase: the huctuations 
of (p are exponentially suppressed and we have (e*"^) 7 ^ 0 at T = 0 or there is a 
quasi-condensation (power-law-decaying order parameter) at T > 0; accord¬ 
ingly magnetic vortex-antivortex pairs become small and dilute, so the gauge 
held is gapped. At the same time the holon, and hence the hole, acquires the 
gap outside the nodes; this is the superconducting phase. Numerically from 
(jlD one can check that for the critical value of A^ the inverse coefficient of the 
Maxwell term in ([2]) is quite small, hence the gauge huctuations are strongly 
suppressed. Therefore the superconducting transition is almost of classical 
3D XY type, driven by condensation of magnetic vortices and related to 
phase coherence as in BEC systems, in spite of the BCS-like charge pairing 
discussed above. This provides a solution of one of the puzzles outlined in 
the introduction. 


3 The loffe-Larkin rule 

We now show a key property in our approach, as well as in the slave-boson 
approach [ 20 ], of the superhuid density of the hole: it can be expressed in 
terms of spinon and holon contributions via lohe-Larkin rule [21], analo¬ 
gously to conductivity. To derive the superhuid density we minimally cou¬ 
ple the holon to an external electromagnetic potential A®™, replacing A by 
A -|- A®™ in the holon contribution of eq. ([2D, denoting the result with 
Ses{p,n^, A^; and we obtain the superhuid density as the coefficient 
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of the quadratic Anderson-Higgs “mass” term in the vector potential in the 
Euclidean effective action for at zero energy-moment a. 

Notice that the space components of fl^ evaluated at 0 frequency-momenta 
give the superfluid density of the holon-gauge subsystem. To obtain from 
n^, Afj_] the effective action in terms of alone we hrst sum over 
the integer currents and integrate over 0. The result can be approximately 
summarized by replacing the XY action of eq. ([3]) by a Gaussian action for 
A, whose coefficient is the superffuid density, p^, of the anisotropic 3D (Vil¬ 
lain) XY model with effective temperature given by 0 = 37r{ms—Al/ms)/Al, 
monotonic in T, with a scale renormalization ^ taking into account the short- 
distance effects. One can now perform the Gaussian integration over A ob¬ 
taining: 


SesiA^P 


1 

2 


d^x 


l /temf rs 




(5) 


The space components of 11^ and those of the kernel of the integral in eq.(|5]) 
evaluated at 0 frequency-momenta give the superffuid density of the holon- 
gauge subsystem and the physical 2D superffuid density ps, respectively. 
Therefore one obtains the loffe-Larkin rule: 


Ps 


PsPs 

Ps+Ps' 


( 6 ) 


A consequence of loffe-Larkin rule in this approach is the 2/3 critical 
exponent. In fact, p^ in BGS approximation is non vanishing below the 
holon pairing temperature Tp^, hence it is still non-vanishing at where 
pI vanishes. Therefore p'l dominates the behaviour of ps near the phase 
transition, and since 0(T) is analytic at Tc, being divergent only at Tps, the 
resulting critical exponent of the superffuid density is that of the 3D XY 
model. 

The theory of the XY model suggests that 

(7) 

and pI{T) = p^(O)pxv(0(7”)/0(Tc)), where pxY is the spin stiffness of the 
anisotropic 3D XY model considered above. In the BGS approximation the 
holon contribution can be derived from the standard formula as 


P^(0) ~ ^ 


d0, , 
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2A^ 
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in PG, taking into account the linear dispersion and the two FS of the 
holons. 
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Figure 1: The normalized superfluid density, for S = 0.12 vs. Our 

theoretical calculation (solid line, ^ = 5.6) is compared with a pure 3D XY 
model (dashed line), data from [7j, £g. 21 (squares), corresponding to the 
near-universal YBCO behavior for superfluid density, from [ 22 ] (circles) for 
5 = 0.075 LSCO and from [23], £g 3b (triangles) for near-optimal-doping 
BSCCO. Data from [ 22 ] and [23] are plotted considering the inflection point 
discussed in the text. The near-universal behavior of our is shown in the 
inset. 

We now compare our results with experimental data. All the parameters 
(except the 5, T-independent scale being hxed by the phase diagram as in 
[9] , we correctly £t the normalized superfluid density data for YBCO (a-axis) 
over a broad range of dopings [211 [T] [251 ES], see £g. [H 

Analogous data in moderate underdoped BSCCO [23], Hg-1201 [22] and 
LSCO [22] are well htted except for a small region near Tc, beyond an inflec¬ 
tion point in the curve, where, however, the authors of [22] claim that doping 
inhomogeneities may affect the reliability of the experimental data. 














4 The Uemura relation 


We now prove analytically that from the temperature dependence of 
follows an approximate Uemura relation. Let us suppose, as the numerics 
proves, that the spinon contribution dominates the superfluid density also at 
T = 0. Since is the solution of a gap equation with critical temperature 
Tps, its temperature dependence will be a function of T/Tpg. Furthermore 
one can check numerically that Tc <C Tps and for a wide range of dopings, 
more precisely slightly away from the doping at which Tps vanishes, A^(T = 
0) ~ Therefore we can parametrize A^(T) = rnlF{T/Tps) for a suitable 
function F obtained from (jl]), which turns out to be almost independent of 
5 and with F(0) ~ 1 in the above quoted range of dopings, see £g. [2], lower 
inset. The effective temperature 0 dehned in the previous section can then 
be written in terms of F as: 


^ Stt 1 - F{T/Tps) 
rus F{T/Tps) 


so that the (dominating) spinon contribution to superfluid density in eq. 
([7]) can also be written as: 


m 


imsTpsFiQ)'^ 

3vr|F'(0)| 


(9) 


If we impose the condition of criticality to the effective temperature of 
the XY model: 0(Tc) = T^'^ 2.2 [27], expanding F linearly, which can 

actually be justihed numerically, we obtain: 


Tp, f (0) - 1 + F(0)mJ/"7(3;r) 

“~|F'(0)| l + m.TXYH3^) ^ 

and the derived shape of Tc{5) in PG is approximately parabolic, as in 
the phenomenological curve hrst proposed in [28], see £g. E) 

In the region where F(0) ~ 1 comparing with ([9]) and using loffe-Larkin 
rule one obtains: 


T 


XY 


-Ps{0) 


(1 - p,(0)/p^(0))(l + msT/^^/{37r)) 


( 11 ) 


Since msT^'^/ (Svr) 1 and the spinon contribution is dominating, eq. 
flTTl) yields the Uemura relation, see £g. E] 
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Figure 2: Theoretically calculated T = 0 superfluid density (solid line, ar¬ 
bitrary units) and critical temperature (dashed line) vs. 6 exhibiting an 
approximate Uemura relation over the broad doping range considered here. 
Insets: approximate (5-universality of the slope of ps near T = 0 (upper inset) 
and of F (lower inset). 


Last but not least, since 0(T) is approximately proportional to T, the 
ratio 0(T)/0(Tc) ~ T/Tc, hence when the holon contribution is small the 
superfluid density ps(T)/ps{0) is approximately the universal function of 
PxviT/Tf.) derived from the 3D classical XY model in the entire temper¬ 
ature range, not only close to the phase transition, and this property is 
quite well verified by the experimental data, as fig. [T] shows. From the 
above property and the Uemura relation one derives also Ps(0) ~ cTc with 
c doping-independent and expanding ps{T)/ps{0) linearly in T/Tc, we can 
write PsiT) — ps{0) ^ c(dpxy/d(T/Tc))(0)T. Hence in the considered doping 
region the T-slope of the superfluid density near T = 0 is almost doping 
independent, see fig. [2], as discussed experimentally in [22] • 
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5 Conclusions 


Summarizing, we have shown that many experimental features of the super¬ 
fluid density of hole-doped cuprates in the intermediate doping region are 
captured by an approach of the hole as a gauge composite of charge and 
spin excitations. With a single {6, T)-independent scale parameter ^ keeping 
track phenomenologically of small scale physics, our model fits rather well 
normalized superfluid density data from YBCO (a-axis), BSCCO, Hg-1201 
and LSCO from moderate underdoping to nearly optimal doping. The uni¬ 
versal critical exponent 2/3 and the near-universality of 3D XY type of the 
normalized superfluid density are also reproduced independently of as well 
as the approximate Uemura relation which is analytically derived. 

These results suggest that the spin-charge composite nature of the hole 
plays a key role in high-Tc superconductivity. 
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